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Abstract. A parameter estimation problem is considered for a one-dimensional 
stochastic wave equation driven by additive space-time Gaussian white noise. The 
estimator is of spectral type and utilizes a finite number of the spatial Fourier 
coefficients of the solution. The asymptotic properties of the estimator are studied 
as the number of the Fourier coefficients increases, while the observation time and 
the noise intensity are fixed. 



1. Introduction 
Consider the stochastic wave equation 

(1.1) — = e, — + e,- + w{t\Q<t<T,Q<x<., 

with zero initial boundary conditions, driven by space-time white noise W . The 
solution of this equation can be written as a Fourier series 

[2 ^ 

u{t,x) = \ — y^Uk{t) sm{kx). 

^ fc=i 

The objective is to construct and investigate the maximum likelihood estimators of 
the unknown numbers > and 02 G M, given {ui{t), . . . , MAr(t)}, t G [0, T], the first 
Fourier coefficients of the solution. 

A similar problem for stochastic parabolic equations is relatively well studied, with 
the first result announced in the paper by Huebner, Khasminskii, and Rozovskii [3]. 
While most of the existing papers concentrate on estimating either a single parameter 
or a function of time, estimation of several parameters in parabolic equations has also 
been studied [21 [8]. The objective of the current paper is to extend some of the results 
from [2] to stochastic hyperbolic equations. In both parabolic and hyperbolic setting, 
the underlying assumption is that the solution u = u{t,x) of the equation can be 
measured at every point in time and space. Then the Fourier coefficients of the 
solution can be computed and used to construct the estimator. 

The main result of the paper is as follows. 
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Theorem 1.1. The (joint) maximum likelihood estimator of the parameters 61,62 is 
strongly consistent and asymptotically normal as N ^ 00. The normalizing matrix is 
diagonal, with the diagonal elements N^^"^ and N^/'^ ; these elements specify the rate 
of convergence of the estimator to 61 and 62, respectively. 

This theorem is proved in Section 3. In Section 2, we estabhsh existence, uniqueness, 
and regularity of the solution of (11.11) . 

Throughout the presentation below, we fix a stochastic basis 

¥ = {Q,J^, {J-ih>o,P) 

with the usual assumptions (completeness of J-'o and right-continuity of J-'t). We 
also assume that F is large enough to support countably many independent standard 
Brownian motions. For a random variable ^, denotes the expectation. M" is an 
n-dimensional Euclidean space; C{A; B) is the space of continuous functions from A 
to B; A/'(m, cr^) is a Gaussian random variable with mean m and variance cx^. 

Finally, for the convenience of the reader, we recall that a cylindrical Brownian motion 
W = W{t), t > 1, over (or on) a Hilbert space if is a linear mapping 

W : f^ Wf{-) 

from H to the space of zero-mean Gaussian processes such that, for every f,g & H 
and t,s > 0, 

(1.2) E{Wf{t)Wg{s)) = mm{t,s){f,g)H. 

li {hk, A; > 1} is an orthonormal basis in H and Wk, k > 1, are independent standard 
Brownian motions, then 

(1.3) f^J2^f,hk)HWk{t) 

k>l 

is a cylindrical Brownian motion. Thus, a cylindrical Brownian motion W is often 
represented by a generalized Fourier series 

(1.4) W{t) = J2Mt)hk, 

k>l 

where Wk = Wh,.- The corresponding space-time white noise is then 

Wit) = J2Mt)hk. 

k>l 



2. Stochastic Wave Equation 



Consider the equation 



(2.1) ^ = a^^-2b^ + W{t),0<t<T,0<x<n, 
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where FT is a cylindrical Brownian motion over L2{{0, tt)). For simplicity, we assume 
(2.2) a > 1, 2\b\ < 1; 

du 



(2.3) u\..„ 



= 0, u\x=o = u\x=Tr = 0; 

t=0 

see Remark 12.31 below about relaxing these assumptions. In physical models, a > 
represents the speed of the wave and b characterizes damping (amplification, if 6 < 0). 

For 7 G M, define the Hilbert space H'^ as the closure of the set of smooth compactly 
supported functions on (0, tt) with respect to the norm 

(2.4) \\fh=(T.^ k^'fi^ > *here h = sj^ f{x) Mkx)dx. 

Note that each of the functions sin(fca;) belongs to every H'^, and if / is twice 
continuously-differentiable on (0, vr) with /(O) = /(tt) = 0, then, after two inte- 
grations by parts, \ fk\ < k^^ sup^(z(^o^T,) \ f"{x)\, so that, in particular, / G H^. More 
generally, every / G H'^ can be identified with a sequence {fk, > 1} of real numbers 
such that J2k>i fk < Even though / is a generahzed function when 7 < 0, 
we will still occasionally write / = f{x), keeping in mind a generalized Fourier series 
representation f{x) = \/2/7r ^^^i fksm{kx). 

Given 7 > 0, / G H^'^ and g G H'^, we define 

if, 9) = ^fk9k] 

k>l 

if G L2((0,7r)), then 

/•TT 

if, 9) = / f{x)g{x)dx. 



In other words, (■, ■) is the duality between H"* and H relative to the inner product 
in = L2((0,7r)); see |5i Section IV.1.10]. 

Equation fl2.ll) is interpreted as a system of two first-order Ito equations 
(2.5) du = vdt, dv = {a^u^^ - 2bv)dt + dW{t). 

More precisely, we have the following definition. 

Definition 2.1. An adapted process m G ^2(^2 x (0,T) x (0,7r)) is called a solution 
of (12.11) if there exists an adapted process v such that 

(1) veL2{n;L,{{0,Ty,H-')); 

(2) For every twice continuously-differentiable on (0, vr) function f = f{x) with 
/(O) = /(tt) = 0, the equalities 

{u{t,-)J)= [\v{tr)J){s)ds, 

(2.6) 

{v{tr),f)= [ {a'{u{tr)J")-2b{v{tr)J))ds + Wf{t) 



4 W. LIU AND S. V. LOTOTSKY 

hold for all t G [0, T] on the same set of probability one. 

Here is the main result about existence and uniqueness of solution of fl2.1l) . 

Theorem 2.2. Under assumptions l \2.2h and (\2.3L equation (12.11) has a unique so- 
lution and, for every 7 < 1/2, 

(2.7) u e L2{Q;L2{{0,Ty,H'')); v G L^iQ; L2{{0,T); H"^-')) . 

Proof. While the result can be derived from the general theory of stochastic hyperbolic 
equations (see, for example, Chow [U Theorem 6.8.4]), we present a different, and a 
more direct, proof. This proof will also help in the construction and analysis of the 
estimators. 



Take in (12. 6p f{x) = y^2/7T sm{kx) and write Uk{t) = (^(t, ■),/), Vk{t) = (f(t, ■),/), 
u!k = Wf. Then 

(2.8) Uk(t) = I Vk{s)ds, Vk(t) = -a^k"^ I Uk{s)ds - 2b Vk{s)ds + Wk{t), 

Jo Jo Jo 

or 

(2.9) ilkit) + 2bukit) + a^k\kit) = Wkit), Mfe(O) = ^^(0) = 0. 
By assumption (12. 2p . 

(2.10) a^k^ > b^ 
for all > 1. Define 

(2.11) 4 = Va^k^ - 62. 

Using the variation of parameters formula for the linear second-order equation with 
constant coefficients, we conclude that the solution of (12.81) is 

Mt) = Y fe-"^'-'^ sin (4(t - s))dwk{s), 

(2.12) 

Vk{t) = -J e'^(*-^' (4 cos (4(t - s)) - 6sin (4(t - s))ywkis). 

By direct computation, there exists a number C = C(T, a, b) such that, for all t, s G 
[0,T], 

(2.13) En^(t) < r,'C{T) = Mit) < C{T). 
Then the Gaussian processes 

(2.14) u(t,x) = \ —''^^Uk{t) sm{kx), v(t,x) = \ —''^^Vkit) sm{kx) 

^ k>l '"^ k>l 

satisfy (12.61) and (12. 7p . Uniqueness of the solution follows from the completeness of 
the system {1/2/7? sm{kx), k > 1} in L2((0,7r)). □ 
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Remark 2.3. We can now comment on the significance of assumptions (12.21) and 
f l2.3l) . Assumption fl2.2p can be relaxed to a > 0, because we will still have a^k"^ > P 
for all sufficiently large k, and so representation formulas (12.121) for the solution of 
equation (12. 9p will continue to hold for all sufficiently large k. In other words, ifa>0, 
then the free motion (any solution of the homogeneous version of (12.91) ) is oscillatory 
for all sufficiently large k > 1; the oscillations are damped ifb > 0, harmonic ifb = 0, 
and amplified if b > 0. This is also the reason to call b the damping coefficient, 
with an understanding that negative damping means amplification. Thus, ( \2.2\\ is 
only needed to simplify the computations by ensuring that equalities (I2.12p hold for 
all k>l. 

Non-zero initial conditions, if sufficiently regular, will not affect existence and regu- 
larity of the solution. Similarly, the analysis will not change much for zero Neumann 
or other homogeneous boundary conditions. 



3. Estimating the Coefficients 

In this section, we assume that the solution u = u{t,x), v = v(t,x) of equation (12.11) 
is observed for all < t < T and x G (0,7r), and study the question of estimating 
the numbers a^,6 from these observations. It will be convenient to introduce the 
notations 

(3.1) Oi = a^, 02 = -2b, 
so that (12. ip becomes 

f)'^ii f^'^ii rill 

(3.2) -W=^'^-^^'m^^^'^' t<0<T, 0<x<vr. 

To simplify the presentation, we keep the assumptions (12. 2p and (12. 3p . By Theorem 
12.21 the solution of (13. 2p has a Fourier series expansion (I2.12p . We will construct 
the maximum likelihood estimators of 6i and 62 using the observations of the 2N- 
dimensional process {ukit), Vkit), k = 1, . . . , N, t E [0, T]} and study the asymptotic 
properties of the estimators in the limit ^ 00. Note that both the amplitude of 
noise and the observation time are fixed. 



By (in, 

(3.3) Uk{t) = I Vk{s)ds, Vk{t) = -Oik!^ / Uk{s)ds + 62 I Vk{s)ds + Wk{t). 



For each k > 1, the processes Uk, Vk, and Wk generate measures P^, P^, P^ in 
the space C((0,T);]R) of continuous, real-valued functions on [0,T]. Since Uk is a 
continuously-differentiable function, the measures P^ and P^ are mutually singular. 
On the other hand, we can write 

(3.4) dvk{t) = Fk{v)dt + dwk, 

where Fk{v) = —Oik'^ Vk{s)ds + 92Vk{t) is a non-anticipating functional of v. Thus, 
the process f is a process of diffusion type in the sense of Liptser and Shiryaev [TJ 
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Definition 4.2.7]. Further analysis shows that the measure is absolutely continuous 
with respect to the measure P^, and 



dP 

(3.5) 



T 

{- 9,k\kit) + 92Vk{t))dvk{t) 



1 '•^ 



{-9ieuk{t) + 92Vk{t)Ydt 



see [TJ Theorem 7.6]. Since the processes Wk are independent for different k, so are the 
processes Vk- Therefore, the measure P'''^ generated in C((0,T);M^) by the vector 
process {vk, k = 1, . . . , N} is absolutely continuous with respect to the measure P'"'^ 
generated in C((0, T); M^) by the vector process {wk, k = 1, . . . , N}, and the density 
is 

■^^(Vk) = exp I ^ y ( - 9ik\kit) + 92Vkit))dvkit) 
(3.6) \fc=i 

( - ^^f^'Mt) + 92Vk{t)fdt ■ 

^ k=l / 



the corresponding log-likelihood ratio is 

ZN{9,,92) = y^[ / {-9,k\k{t) + 92Vk{t))dvk{t) 



(3-7) 

1 '•^ 



{-9ieuk{t) + 92Vk{t)Ydt 



Introduce the following notations: 

N „T N „T 



(3.8) 



Ji,N = Y^k^ ul{t)dt, J2,iv = X^/ ^lit)dt 

k=l k=l 
N T 

Ji2,N = y2^^ '^k{t)vk{t)dt; 

k=l 

N T ^ 

Bi,N = - ^k^ / Uk{t)dvk{t), ^i,N = ^k'^ / Uk{t)dwk{t); 

k=l "^0 k=l "^0 

B2,N = X] / Vk{t)dVk{t), ^2.N = / 
fc=l "^0 fc=l 



Note that the numbers J and B are computable from the observations of Uk and f^, 



A; = 1, . . . , A^, and also 



(3.9) Bi^N — 9iJi^N — 92Jl2,N — ^l,N, B2,N — —9lJl2,N + 92J2,N + C,2. 

1 ^ 

(3.10) Ji2,iv = :T5^A;2n^(T). 



N, 



2 

fc=i 
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We consider the problem of estimating simultaneously both 9i and 6*2 from the ob- 
servations 

K(t), Vk{t), k = l,...,N, te[o,T]. 
The maximum likelihood estimators 6i n, ^2 v satisfy 



and 



dZ]\f{6i, 62 



0, 



or, after solving the system of equations, 

in h BinJ2^N + B2^nJi2,N ^ 

{6. II) Ui^N — J J 72 ' '^S.Af 

Ji,nJ2,N — ^12 N 



BinJi2,N + B2NJ1, 



N 



Ji,nJ2,N — J12 N 



For T > and 02 E M, define 



(3.12) 



Cie2,T) 



.02T _e^T-i 



if 02 + 0; 
if ^2 = 0. 



Note that C(02, > for all T > and ^2 e M. 

The following theorem describes the asymptotic behavior of the estimators fl3.11l) . 
Theorem 3.1. Under assumptions fl2.2l) and fl2.3p . 



hm 9^^], = 9,, hm 9^^^ = 92 



with probability one and 



39i 



in distribution. 



Proof. Define 



D 



N 



2,N 



It follows from ([S3D and fIXTTD that 



(3.13) 



^2^Af — O2 + 



1 



1 — Dn \Jl,N 
1 / ^2,N 



^1,N . ^ Jl2,N 
+ 'i2,N 



+ 6, 



Ji,nJ2,N 
Jl2,N 



N- 



^ — B>N \ J2,N ' Ji,nJ2,N , 

By direct computations using fl2.12p (and keeping in mind (13. ip ). 
(3.14) 



hm rul{t)dt= ^^^''^\ 
Jo Oi 



k—>oo 
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and 

(3.15) lim iV-3EJi,;v = 

Since each Uk is a Gaussian process, 



sup fc^E / ul(t) dt < oo, 
k Jo 



and then the strong law of large numbers implies 

lim — - — = 1, hm 



N-*oo EJi^AT N^(x> EJi^Tv 

both with probability one [apply the first theorem in Appendix, taking = k'^ Jq u\dt 
and then, = k'^ Uk{t)dwk{t)]. The central limit theorem implies 

lim ^ ^''^ = MO, 1) 

in distribution [apply the second theorem in Appendix, taking fk{t) = k'^Uk{t)]. Sim- 
ilarly, 

(3.16) limE/ vl{t)dt = C{92,T), 

Jo 

and 

(3.17) lim N-^EJ2,N = 0(02, T). 
Since each Vk is a Gaussian process, 

supE / vl{t)dt < oo, 
k Jo 

and then the strong law of large numbers implies 

lim ^'^ = 1, lim ^^'^ = 0, 

N->oo EJ2,Ar N^oo EJ2,7V 

both with probability one. The central limit theorem implies 

= Ar(o, i; 

in distribution. Finally, define 



I if ^2 ^ 0; 

cie2,T) = { ^202 

r, if 02 = 0. 



Then ( I330D and ([212]) imply 



hm Ar-iEJ,2,^ = 
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and, by the strong law of large numbers, 

,. Jl2,N 

lim = 1 

N^oo EJi2,Ar 

with probability one. Then fl3.15p and fl3.17p imply 

Jl2,N C{(^2-,T) 



lim = 0, lim 



JV-oo " ' N^oo J2,j^ 2610(62, T)' 

both with probability one. The conclusions of the theorem now follow. □ 
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Appendix 

Below, we formulate the strong law of large numbers and the central limit theorem 
used in the proof of Theorem I3.1[ 

Theorem 4.1 (Strong Law of Large Numbers). Let C,^, k > 1, be independent random 
variables with the following properties: 

• = 0, > 0, 

• There exist real numbers c > and a > —1 such that 

lim k-'^E^l = c. 

k—*oo 
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Then, with probability one, 

lim %^=^ = 0. 

If, in addition, E^^ < Ci ^E^l j for all k > 1, with Ci > independent of k, then, 
also with probability one, 

lim 2^^^!^ = 1. 

Proof. This is a particular case of Kolmogorov's strong law of large numbers; see, for 
example, Shiryaev [9, Theorem IV. 3. 2]. □ 

Theorem 4.2 (Central Limit Theorem). Let Wk = Wk{t) be independent standard 
Brownian motions and let fk = fkif) be adapted, continuous, square-integrable pro- 
cesses such that 

in probability. Then 

lim ^^-^^°"^^^^^^"^^l =Ar(0,l) 



xti^io mdt 

in distribution. 



Proof. This is a particular case of a martingale limit theorem; see, for example Jacod 
and Shiryaev [4, Theorem VIII. 4. 17] or Liptser and Shiryaev [6", Theorem 5.5.4(11)]. 

□ 
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